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ABSTRACT 


Ceramic matnx composites are promising materials for applications in a wide 
temperatuie range, fiom ambient to a very high tempeiature Many of such 
composites aie known foi then excellent mechanical and chemical properties at an 
elevated temperature In addition, these matenals are light weight and retain their 
mechanical strength at high temperatures 

In the recent past, chemical vapor infiltration (CVI) has emerged as a 
piomismg method for processing ceramic matnx composites, in which the matnx is 
chemical vapoi deposited m a fibrous prefoim 

Piocessmg of ceramic matrix composites depends upon many piocess 
vaiiables, such as tempeiature, temperature gradient, pressure, flow rate of piecursoi 
gas, structure of pieform (that is, natuie of poie netwoik), and concenti ation of 
piecuisoi gas 

In the present work a model is developed foi diffusion m porous media and it 
is applied to study the effect of pore orientation m the deposition of SiC preform by 
SiC matrix 

The model uses the Dusty gas model tor descnbing diffusion which takes into 
account binary diffusion, Knudsen flow, and viscous flow The concept of onentation 
averaging foi diffusion of species in a poious media is applied The poious medium 
is divided into many elemental volume each containing pores oriented m a specified 
direction and by applying the orientation averaging method the porous medium is 
repiesented by a piece wise homogeneous medium 

We considei foui cases of pores orientations at angles of <j> = 0°, 45°, 90° and 
135° and the various combinations included pores in single, two, three and four 


in 



directions The studies are carried out at a temperature of 1200 K for the initial pore 
radius of 2 5 |im The initial porosity for different pore orientations ranged from 0 3 
to 0 45 

The results show that open-ended pore and pores inclined to the surface of 
the preloim enhance the densilication of the preform 
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CHAPTER 1 


INTRODUCTION 


Ceramic matrix composites are promising material for applications in a wide 
temperature range, from ambient to a very high temperature Many of such 
composites are known for their excellent mechanical and chemical properties at an 
elevated temperature In addition, these materials are light weight and retain their 
mechanical stiength at high temperatures Thus, they find applications in heat 
generators, turbines and space systems[l] 

Processing of ceramic matrix composites can be classified into two categories 
In the first category, the fiber-matrix preform is prepared by slurry infiltration or sol- 
gel method and then the preform is either hot pressed or sintered In the other 
category, the composites are prepared by melt infiltration or chemical vapor 
mfiltiation (CVI) 

In the slurry infiltration method, the fiber tow is passed through a slurry 
mixture, constituting of a carrier liquid, matrix powder and a binder, which is usually 
an organic compound The tow is subsequently dried, cut and stacked one above the 
other in a desired configuration Then, the composite is made by sintering or hot 
pressing the stack In the sol-gel process, the fiber is dispersed in the sol The gel 
material is allowed to dry to form a porous ceramic, which is then densified by 
sintering Both the methods suffer from serious problems A major disadvantage of 
the slurry process is damage to the fiber during hot pressing, where the temperature is 
high On the other hand, sol-gel processing involves a large volume shrinkage during 
drying and also during densification firing stage Furthermore, the yield of ceramic 


material is low 



These disadvantages can be overcome by processing ceramic composites by 
melt infiltration or chemical vapor infiltration However, melt infiltration is of a 
limited utility in producing ceramic matrix composites, because it requires melting of 
the matrix material, which in most cases melts at a very high temperature 

In the recent past, chemical vapor infiltration (CVI) has emerged as a 
promising method for processing ceramic matrix composites In this process, the 
matrix is chemical vapor deposited in a fibrous preform Thus, the process is a 
specialized form of chemical vapour deposition (CVD) process, except that here, 
instead of gas flow over the surface, the precuisor gases, diffuse into a porous preform 
made up of fibers Chemical reaction at the pore wall leads to the deposition of the 
matrix material 

The chemical vapor infiltration process has many advantages over the other 
processes Fabricating composites using this process requires only atmospheric 
pressure and low temperatures compared to other processes, thus minimizing damage 
to the fibrous piefoim 

1 1 LITERATURE REVIEW 

Fiona heie on we focus on fabrication of SiC fibei-SiC matrix composites by 
CVI This composite has a high strength and elastic modulus, and exhibits good 
mechanical strength associated with high strain-to-failure Also, SiC composites are 
resistant to high temperatures, and have high thermal conductivity, and low thermal 
expansion, which provides good thermal shock resistance^, 3] 

Trichloromethylsilane (TMS) is used as aprecursoi gas in fabricating 
S 1 C/S 1 C composites[3, 4] The matrix is deposited by thermal decomposition of TMS 
in excess of hydrogen [5] The decomposition reaction proceeds according to 



CH3S1CI1 (g) = SiC (s) + 3 HC 1 (g) (in excess of H2) 

Chemical vapor infiltration process for fabricating S 1 C/S 1 C composites can be 
classified into four categones[l] 

Isothermal chemical vapor infiltration (ICVI)[4, 6] 

Thermal gradient chemical vapor infiltration [7] 

Forced flow chemical vapor infiltration (FCVI)[8, 9] 

Pulse flow chemical vapor infiltration[10] 

In isotheimal chemical vapor infiltration (ICVI), the preform is maintained at 
a uniform temperature in a chamber and the precursor gas is flown over the 
preform[4] The major disadvantage in this process is that, at high temperatures, a 
preferential deposition occurs on the surface layer, leaving behind a large internal 
porosity In order to reduce the residual porosity, if the process is conducted at a 
lower temperature, it takes longer time to densify the preform 

This difficulty is overcome in thermal gradient CVI, where a temperature 
gradient is applied across the thickness of the preform, and the gases diffuse through 
the low tempeiature end[7] This procedure allows densification such that deposition 
front move from high temperature end to the low temperature end 

In FCVI[8], the precursor gas is forced through the preform in which a 
temperature gradient is also maintained In yet another process, as in the previous 
case called pulse flow CVI, pulsed flow of precursor gas is caused by cyclic 
evacuation of the reaction chamber and back filling with precursor gas [10] 

Despite the mentioned disadvantages for ICVI, it is frequently used because a 
large number of parts of different sizes and shapes can be easily made and also very 
irregular shapes can be fabricated using this process This process is commercially 
used in the fabricating high performance brakes for aerospace applications^] 



In the following, I review various models for diffusion in porous media and 
chemical vapor infiltiation process related to SiC composites The models for 
diffusion of species in porous media can be broadly classified into two approaches, 
discrete and continuum In the first approach, the porous media is treated as a discrete 
collection of objects, such as voids and grains In the continuum models, porous 
media is treated as a continuum were temperature, concentration, and solid species 
concentration vary smoothly with respect to position and time[l 1] 

Among the discrete approaches, use of percolation theory is common Reyes 
and Jenson [12] simulated pore related properties in porous media using the concept of 
percolation theory on a bethe lattice They evaluated effective transport coefficients 
without introducing a toitousity factor Their model inherently accounted for narrow 
necks, tortuous path and dead ends 

Another approach called “Effective medium approximation (EMA),” was 
originally applied by Landauer[13] to find the effective electrical conductivity of a 
bimetal-strip Subsequently, Kirckpatrick[14] used this theory to find electrical 
conductivity of a random resistor network 

Burganos and Stonchos[16] applied EMA (or EMT) to find the effective 
conductance of the pore network and, by assuming that network of pores satisfies the 
smooth field approximation (SFA), they applied these conductivities to find effective 
diffusivity in the pore network Sotirchos[17] further developed the model by adding 
convection along with diffusion in a porous media Pedro et al [18] then extended the 
EMA theory for anisotropic lattice structures 

Two main appioaches used in continuum models are dusty gas model and 
capillary model In “dusty gas model”[ 19-21], the solid is treated as another heavy 
immobile gas, which is uniformly distributed in the medium By treating the medium 



as a heavy immobile gas, transport properties are calculated by applying the kinetic 
theory to all the gas species including the “dust gas” The binary diffusion coefficient 
between the dust gas and another gas particles is equivalent to Knudsen diffusion 
coefficient Further the viscous flow arising from pressure build up is also included in 
the mass flux expression Thus, the flux expression includes all, the multicomponent 
ordinary diffusion, Knudsen diffusion and viscous flow Therefore, the model applies 
well to diffusion in the transition regime However, the effect of structure of the 
porous medium is included only as empirical parameters in molecular, Knudsen and 
viscous flow 

Though dusty gas model gives satisfactory results, its ability to deal with the 
structure of porous medium is not fully developed 

Another important model in the continuum approach is Feng and 
Stewart’ s[2 1-22] capillary model They represented the medium by a network of 
pore, which is thoroughly cross linked By assuming the species concentration and 
pressure to vary smoothly with respect to position and time, they evaluated the flux of 
a species for a fraction of all pores aligned in a particular orientation and of a given 
radius Then, the total flux for a species was calculated by integrating the flux over 
all orientations and radii However, the flux expression was based on dusty gas model 
and, in addition, the model also included surface diffusion in the flux expression 
Thus, Feng and Stewart’s model was the first attempt to include the effect of pore 
network in the dusty gas model 

Apart from capillary model, grain model is also used where solid phase 
consists of sphencal paiticles or grams Evangelos et al [23] modeled a gas-solid 
reaction with partially overlapping grains, each consisting of a hard core and a soft 
shell in which fluid could penetrate or diffuse Only the soft shells were allowed to 



oveilap with each other as the reaction proceeded This approach was verified their 
mode] by comparing with the experimental data for the reaction of ZnO particles with 
0 5 % H 2 S in N 2 at 600°C 

Another model involves volume averaging of transport properties over a 
domain involving solids and voids Using the volume-averaging model Jin et al [24] 
studied the phenomenon of diffusion in the media consisting of glass sphere, mica 
particle and disks made of mylar sheets 

Another fiequently used model to describe porous media is pore network 
model Hollewand and Gladden[25] modeled a random network of pore to study the 
problem of diffusion and reaction occurring simultaneously in a porous catalyst 
Tortuosity is a common parameter in modeling a porous media and result from this 
model suggests that the tortuosity obtained under strong diffusion control cannot be 
applied to the condition where both diffusion control as well as reaction control is 
important In another work the same authors compared random and regular network 
using percolation and finite scaling theory[26] and found that random pore model had 
gi eater percolation threshold They also developed a random micro and macro pore 
network by connecting the micropore to macropore and obtained experimentally 
matched toituosity Reckmann and Keil[27] modeled a three dimensional pore 
network with piedefined poie radii, connectivity and porosity for any type of reaction 
Dusty gas model was used for calculating mass transport in a pore Fuertes and 
Marben[28] used random pore model, incorporating percolation phenomenon to 
analyze gasification reaction when percolation takes place simultaneously Zhang and 
Seaton [29] compaied the diffusion coefficient obtained by continuum model to that 
from pore network model to investigate the effectiveness of using continuum 
diffusion -leaction model to descnbe diffusion and reaction in pore network When 



the pore network was very close to percolation threshold, the continuum diffusion- 
reaction model failed 

Pavol and Yang[30] used bethe lattice to represent the porous media They 
used percolation model of mass transport in the network and applied effective 
medium theory to find the effective diffusivity of condensable vapor in porous media, 
where condensation and adsorption-desorption hystenes occurred 

In addition to various models to treat transport in a porous media, we now 
provide a review on modeling of chemical vapor infiltration of SiC Gupte and 
Tsamapoulos[31] modeled mass transport and deposition of SiC by chemical vapor 
infiltration, in a single pore Their model predicts negligible effect of pressure on 
deposition as long as the pressure is below 0 5 atm and a temperature less than 
1000°C, if isothermal chemical vapor infiltration is used The authors also modeled 
deposition using temperature gradient and found that small temperature gradient 
enhances the densification piocess 

Lu[32] studied the effect of a structural parameter and theile module on 
densification using a random pore structure model In that work, the structural 
parameter is a non-dimensional quantity which depends on surface area of the pore, 
porosity and total poie length He found that a smaller value of theile module leads to 
moie unifoim deposition, whereas pore plugging occurrs at high theile module and 
stiuctural parametei Furthei, substance having wide pore distribution was found to 
offer more resistance to diffusion as deposition proceeds (that is after initial 
deposition) 

Chung and smith[33] modeled deposition of SiC by chemical vapor infiltration 
on a sample of overlapping fibrous preform Their model assumed diffusion to occur 
in channels formed between the fibers, space between stacked mats and channel 



formed due to stacking of mats one on top of another Comparing two cases in model, 
where in one channels are aligned straight and are staggered in the other, the former 
alignement leads to a lower porosity 

Deepak and Evans[34] modeled deposition of SiC from tnchloromethylsilane 
in a porous preform in which thermal gradient was achieved by microwave heating 
Dusty gas model was used for describing mass transport in the preform They suggest 
that at least one dimension of the preform should be short for gas to diffuse through 
the preform Lower temperature, and larger pore size favored uniform deposition of 
matux material, but required longer time for densification Also temperature gradient 
favors more densification in less time Kawasi[35] et al fabricated SiC whisker 
reinforced alumina composite by chemical vapor infiltration technique, where the 
core was heated to 1 100°C at a total pressure of 3 3-13 kpa Preform surface was 
cooled by radiation as well as convection due to gas flowing around the preform 
surface The authois also modeled isothermal and thermal gradient CVI by 
employing cubic airay of whiskeis Pore plugging occurred in both the cases, but 
isotheimal process had premature pore plugging compared to thermal gradient 
process Comparing both the cases, thermal gradient process greatly reduced the 
processing time 

Chang et al [36] studied the effect of reaction mechanism in the chemical 
vapor infiltration for production of SiC composites by consideung different reaction 
model for deposition of SiC They inferred that first-order kinetic model is not 
correct for predicting the infiltration process by showing influence of chlorine 
containing by-pioducts in delaying the infiltration time Use of alternate method of 
precursor chemistry, which avoids chlorine containing by-products, was suggested for 
minimizing infiltration time Tai et al [37] also considered the influence of deposition 



mechanism in FCVI Two different reaction mechanisms were considered for 
deposition A Zero-order reaction mechanism with activation energy of 188 kJ/mol 
and a first —order reaction mechanism with activation enrergy of 120 kJ/mol were 
adopted for comparison The deposition mechanism was found to effect the density 
distribution in the preform More uniform densification was observed when first - 
order reaction mechanism was considered, but rapid densification occurred for the 
case of zero-oider mechanism They concluded that assumption of reaction 
mechanism affects not only the infiltration time, but also the theoretical matrix 
deposition 

Ofon and Sotnchos[38] studied the effect of multicomponent mass transport 
on the densification of S 1 C/S 1 C composites, using dusty gas model for describing 
transport in the fibious preform They compared various variation of dusty gas model 
by neglecting the intei action of TMS, HC1, and H 2 for two cases, one by including 
viscous flow and the olhei by neglecting the viscous flow The study showed that 
interaction of HC1 with TMS influences the deposition, as diffusion of TMS into the 
preform is impeded by diffusion of HC1 out of the preform It was also demonstrated 
that full dusty gas model should be used unless both temperature and pressure are 
low The same authois also simulated 2-D, 3-D CVI models for densification of 
preform with uniform porosity and an amsotiopic preform (Fibers oriented in one or 
two dimensions) for studying the effect of preform anisotropy, aspect ratio and shape 
of the preform [39] They observed difference in deposition for anisotropic structure 
(having cubic shape) m different direction when it reached percolation threshold For 
large aspect ratio of the preform, deposition in the preform was influenced by gases 
diffusing through the small sides of the preform, whereas the gases diffusing along the 
longer length of the preform effected only the outer surface Shape of the preform 



also influenced the deposition of the gas in the preform, more uniform deposition was 
seen in regular shape(Like cube) Thus they favored multidimensional model to be 
used for simulation unless mass transport is identified to be occurring in a dominant 
direction 

The effect of pore structure during densification of S 1 C/S 1 C composites was 
also studied by Ofon and Sotirchos[40] Two models were considered, one in which 
poies were represented as cylindrical capillaries randomly oriented in all the three 
dnections and another the solid phase was represented by fibers oriented in one, two 
and three dimensions In the later case fiber orientation and direction of diffusion 
were found to effects the overall diffusion coefficients and percolation threshold(that 
is minimum porosity left which cannot be filled) For the case were capillary were 
used to model the porous media, the deposition rate was dependent on the pore size 
distribution The model suggested that use of bi-modal capillary structure (that is, 
two different sizes of capillaries were used to represent porous media) is more 
applicable for simulating a chemical vapor infiltration process The effect of mass 
transport on optimal pressure and temperature pair that can be used for minimum 
infiltration time was also studied by the same authors[41] The optimal pressure for 
CVI under condition of diffusion driven transport was found to be larger than that 
used in the industrial CVI piocess Optimal pressure was found to be located at the 
molecular diffusion regime, and as a result, it tends to increase with decreasing pore 
size It was found that optimal temperature for densification decrease significantly 
with increasing deposition uniformity in the final product, final conversion level at the 
external suiface, reveise rate constant and the preform size Mole fraction of the 
piecursor gas also was found to have a strong influence on processing time 



1 2 OBJECTIVE OF THE PRESENT STUDY 

Processing of ceramic matrix composites depends upon many process 
variables, such as temperature, temperature gradient, pressure, flow rate of precursor 
gas, structure of preform (that is, nature of pore network), concentration of precursor 
gas Optimizing these process variables not only reduces the infiltration time, but also 
affects the structure of matrix material deposited, ultimately affecting the properties of 
the composite 

Piefoim for fabrication of composite material consists of pores, which are to 
be filled by CVI Although at microscopic level, the diffusion is isotropic, but at this 
level it is not possible to model diffusion and deposition because this scale is shorter 
than the pore radius Therefore, models are built at a scale much larger than the pore 
radius However, then at this scale, the effect of pore network make diffusion 
apparently anisotropic 

In the piesent work, I propose to model a porous medium by representing the 
pore network as a piece wise homogeneous medium The basic idea is to use 
orientation averaging method, similar to that by Feng and Stewart[22] But this 
model differs from the model of Feng and Stewait in the way the porous media has 
been represented Feng and stewart represented the whole domain as a homogeneous 
medium, wheieas I will be dividing the porous medium into a number of smaller units 
and by performing onentation averaging in each volume element I will represent the 
pore network as a piece wise homogeneous medium As an application of the model 
developed, I apply the model to a S 1 C/S 1 C system 
The objective of the present study is to 



Repiesent a porous media consisting of network of pores aligned in various 
orientation by an effective homogeneous media consisting of finite number of 
elemental volume 

Study the effect of orientation of network of pores on the deposition of SiC in 
a SiC preform by chemical vapor infiltration 



Represent a porous media consisting ot network of pores aligned in various 
orientation by an effective homogeneous media consisting of finite number of 
elemental volume 

Study the effect of orientation of network of pores on the deposition of SiC in 
a SiC preform by chemical vapor infiltration 



CHAPTER 2 


MODEL FOR CHEMICAL VAPOR INFILTRATION FOR 

SiC COMPOSITES 


This chapter discusses mass transport and chemical reaction in a porous media 
constituting of a capillary network by orientation averaging method The method is 
used with a goal to represent a porous medium with an effective homogenous 
medium In this model orientation averaging is done on a second phase (pores m 
poious media) and an effective equivalent homogenous medium is represented 
through an anisotropy tensor developed in the following sections 

2 1 DIFFUSION IN A HOMOGENOUS MEDIA 

Consider diffusion of v gaseous components in a homogenous volume /? If 
we define diffusion velocity of species l in an arbitrary coordinate system as, 

( 1 ) 

V, = v, - v 0 


The expiession foi V„ based on Chapman-Enskog theoiy of gases [42] is 
derived to be, 


, 2 V 


v, = — 

n tP H 


(2) 


in which thermal diffusion is neglected and where the quantity dj is defined as, 



dj =V Xj + 


S n j*V 


„ n ,m, 
VlnP- J J 


Pp 
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Fj - ^n k F k 


(3) 


Further, we can write the mass diffusion flux from equation (2) as, 


n2m , v ^ „ 

J. = — —2>A d j 


(4) 


Coupled with an equation for total momentum balance (such as the Navier- 
Strokes equation) and an appropriate boundary condition on dj3 , the fluxes can be 
determined in the volume /? However, the problem is much more complicated when 
the volume /? consists of pores rather than a homogenous medium because the 
diffusion occuis preferentially along the direction of the pores But, it is still possible 
to solve equation (4) if the structure of the porous media were known in complete 
detail Unfortunately, it is not always practical to solve the fluxes on the same 
dimensional scale as the pores 

Therefoie, the purpose of the following development is to piece wise 
homogenize the porous media occupying volume, in order to account for the 
amsotiopic natuie of the diffusion caused by the structure of the porous media 


2 2 MATHEMATICAL DEVELOPMENT 

Consider a volume /? which is riddled with cylindrical pores, where a pore is 
defined as a void with one dimension much longer than the other two We subdivide 
the volume into N smaller regions /? a , such that, 



and the radii of the pores in each volume element is same, while in order to take 
advantage of homogenization, variation in the pore radius is allowed from one 
element to another 

We define two coordinate systems [x] and [x / ], related to each other through 
rotation of axes by an angle 0 and (f) as described in Appendix (A) [x] is the global 
coordinate system in which the diffusion flux is expressed after homogenization 
(averaging) Now, consider a pore in with its longest dimension (length) along 
[x 3 ] Thus, [x'] is the local coordinate system associated with a pores and, hence, 
molecular diffusion flux given by equation (4) can be easily expressed in this 
coordinate system Since the other two dimensions (x, andx' 2 ) of the pores are short, 
it is reasonable to neglect the gradients of mole fractions, and pressure in these 
directions 

Thus, equation (4), the diffusion flux in a pore is expressed as, 
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Intioducing L allows foi a non-zero flux only in the [x' 3 ] direction 



The primed quantities are expressed in [x'] and the superscript p on the 
diffusion flux is used to designate that this flux is only in one pore 

We now transform the vector quantities expressed in the local coordinate 
system of the pore, [x ] , to the global coordinate system [x] Thus, according to the 
transformation, equation ( 5 ) is written as 
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( 7 ) 


and dj is now expressed in [x], and the mass diffusion flux, jf , represents the 
diffusion of a component 1 in a pore oriented in the (0, <jj) direction (subsequently 
called direction g) In the following section, a piocedure based on orientation 
averaging is developed to add the contribution, weighted over an orientation 
distribution function, from many pores such that /? a can be treated as an effective 
homogenous media 

We now represent the mass flux in the Cartesian coordinate system Thus, it is 
reasonable to take / 3 a as a cuboid with Lj, Lo and L3 as the length of the three sides 
along the x, ,x 2 ,andx-, direction and Ai, A2 and A3 as the area of the faces 
perpendicular to these three axes, respectively Furthermore, suppose there are 
dN(g, i) pores in orientation [g, g + dg] and of length (/, i + df ] Then, the rate of 
mass diffusion is,j 1 p (g,^)adN(g,f), where ‘a’ is the cross sectional area of each pore 



Thus, the mass diffusion flux, expressed over the superficial area of a volume 
element, is given by 


f 1 
L] A, 

dj 1 (g,^)=adN(g,£)^ 0 

0 


L, 

l 


0 

1 l 2 

L 2 A 2 t 

0 


0 


0 j, p (g,/) 

1 l 3 
l,a 3 l . 


( 8 ) 


Where, l is the length of the pores in the orientation g Also, since LiAi = L 2 A 2 = 
L^A-i is the volume of /3 a equation (8) is therefore, rewritten as, 


dj, (g 5 ^)=de(g,^)Pj 1 p (g, 


(9) 


where, 



is a dyad which accounts for reduced (increased) diffusion due to longer (shorter) 
pores 

In writing this equation, we have assumed that a pore both begin and ends at dj3 a 
Thus only those poie are considered which are connected 



2 3 ORIENTATION AVERAGING 


We define an onentation distribution function such that /( g,^)dgdl? is the 

fraction of porosity in dg onented along g and of length between l and f +d£ 
Expressing it mathematically, we have 

oo) 

£ 4n 

and since the distribution function is normalized, therefore, it satisfies the following 
equation 


4% 


2hk £ 


JJJ/CM , ^)sin 0 d^d9 d(j) 


: 1 


(H) 


0 0 0 


Thus, using equation (9), we write the total flux in /?« expressed over its 
superficial area as, 


2kk t 


(12) 


J,=sjj}Pj , p (e^)/(0,4),^)sm0dMed(j) 


0 0 0 


In order to evaluate the integral in equation (12), we make an assumption that 
in /?a the term, dj appearing in equation (7) is constant Thus, an upper limit on the 
size of is detei mined by ensuring the validity of this assumption However, it must 
be noted that the flux is calculated this way only in a /? a and that these fields may vary 
from one volume element to another Further, after combining equation (7) and (12), 
the mass diffusion flux in the homogenized volume is given by 
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(13) 



where, 


£ 2 KK * 

(a) =— J [Jp (M T LM)f(9,(j)j)sinQ dfdBdcj) ^ 

™ 000 

is the anisotropy tensor 

The anisotropy tensor, which is expanded in appendix (B) and the mass 
diffusion flux in an effective homogenous volume, derived in equation (13), applies to 
any /? a Hence equation (13) and (14) represent the transformation of mass flux m a 
porous media to piece-wise homogenous media of dimension L|, L2 and L3 along the 
thiee respective axes An advantage of this procedure is obtained when the scales Li, 
L 2 and L-) are much larger than the spacing between the pores Yet a J3 a must be kept 
small to the extent that the concentration gradient remain constant within the /?« 

The multicomponent diffusion coefficients in equation (13) are a spatially 
dependent quantities Moreover, the value of diffusivity reported in literature is 
generally for the binary diffusion coefficient Although the two diffusion coefficients 
are related, because of the reasons, given above, it is often inconvenient to use 
equation (13) 

Hirschfelder et al [42] have shown that the equation (4), even in the presence 
of non-isobanc diffusion can be manipulated to obtain a set of generalized stefan- 
maxwell equation, which contains the binary diffusion coefficient instead of the 
multicomponent diffusion 

Similarly in presence of ^A^, equation (13) becomes,- 


iy — — 

n j irqnij'O, 


■(x I m I j J -Xjmjj, 



(15) 



The equation above represents (v - 1) independent equation in the mass fluxes 
j and describe multicomponent molecular diffusion in a porous media However, 
when the pore radius is small, in addition to the molecular diffusion, the viscous flow 
originating from a pressure buildup and the Knudsen flow are important A model 
that includes mass transport by these mechanisms is “dusty gas model” An equation 
similar to the dusty gas model can also be obtained after orientation averaging 
Following the same method as observed by Mason and Malinauskas[20], from 
equation (15), we obtain 
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Also, after using a piocedure similar to onentationally averaging the diffusion 
flux, the viscous flux, based of fully developed Poiseuille flow, is given as 
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(17) 


2 4 APPLICATION TO CHEMICAL VAPOR INFILTRATION 

This section deals with two-dimensional modeling of chemical vapoi 
infiltration (CVI) for production of S 1 C/S 1 C composites Orientational averaging 
method discussed in the previous chapter with “dusty gas model” is used to expiess 
the molai fluxes 

First a porous prefoim made up of largely SiC fibers is assumed Then the 
preform is densified by CVI, wheie in, the gaseous reactant, 
lnchcloiomethylsilane(TMS), and the cainei gas, H? diffuse into the pores 
Following chemical decomposition of TMS according to equation (18), 



CH 3 SiCl 3 (g) = SiC (s) + 3HC1 (g) 


(in excess Ho) 


(IB) 


SiC is deposited on the pore walls The product gas HC1, diffuses out of the pores 
In Fig 2 1 , we show that a two dimensional region is divided into many 
volume (area) elements Ax of dimension Li and L 2 in x and y direction, respectively 
Because we model two dimensional arrangement of pores, the third dimension, L 3 , 
could be of any magnitude /3 a contains a network of pores of specified orientation, 
and volume of pores over which chemical reaction takes place is related to volume of 
Ax (V a = LixL 2 xL 3 ) as 

V a p (t)=s a (t)V a (19) 

And the surface area over which the reaction takes place is given by, 

- ' r| V 

In writing equation (20), we have assumed purely cylindrical pores, neglecting 


the effect of pore intersections 




2L X 


Fig 2 1 Schematic diagram of a preform divided into 
elemental volume 


As the deposition in SiC is described by fust-order, heterogeneous kinetics, 
the deposition of SiC on poie walls in a /3 a occurs according to, 


_i_^c =kcX| , (21) 

~ s « dt 


where k = k 0 e 



and the subscript 1 stands for TMS 


Making use of equation (4) in a mass balance on TMS we get 


N‘ y L|L, +N[^L 2 L 3 — N,“ L ( L, -Nf x L 2 L 3 =S a kcX, p 


( 22 ) 


N 



Fig 2 2 Schematic sketch of a control volume 

Where Lj, L 2 and L3 are length, width and the height of the control volume 
respectively and Ni is molar flux of TMS The supei script indicates the location at 
which flux is considered (see Fig 2 2) The equation (22) is solved using the finite 
difference method 

The multicomponent diffusion flux at the points such as p is given according 
to the dusty gas model, which includes ordinary diffusion, viscous flow and Knudsen 
flow Detail of this flux expression is available in monograph on the Dusty gas 


model [20] 



In the present analysis, the contribution of thermal diffusion is considered 
negligible In that case, the dusty gas model based expression for flux for the 
effective medium in equation (16) is converted to molar fluxes and that reads 
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( 23 ) 


for i,j =1,2,3 

where the subscript 1, 2 and 3 is reserved for TMS, HC1 and H 2 , respectively, and D,* 
is the Knudsen diffusion coefficient for species 1 The quanities i3,j, D,k and r\ are 
estimated by the formulae given in Appendix (C) 

With an assumption of an ideal gas, the following equation of state holds 
P= cR g T (24) 

Fuithermoie, from chemical reaction in eq (18), for every one mole of TMS 
consumed, 3 moles of HC1 is produced Therefore in quasi-steady state, we have 
N 2 = -3N, (25) 

As H 2 does not take part in the reaction , H 2 is stationary and hence 
= 0 (26) 


After explicitly writing eq (22) for species 1, 2 and 3 and then substituting 
eqs (25) and (26), we get three equations, which are 
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( 28 ) 




X,+X 2 +X,-l (30) 

The eqs (22), (24), and (27)-(30) constituting a system of partial differential 
equations in vanableX, ,X 2 ,X 3 ,Pandc are solved subject to the boundary conditions 
discussed below 


For a preform symmetrical with respect to x and y axis, calculation are 
reduced to only the first quadiant and the boundary conditions are, 
at x = L x or y = L y 


X,=Xf, 

x,= x°, 

X,=x°, 

c=c°, 

P = P° 


Also, the boundary condition at x=0 is N v = 0, and similaraly, the boundary condition 
at y = 0, is N x = 0 

When the entire preform is considered due to lack of symmetry, the boundary 
conditions are, 

At x= -L x and L x or y = -Ly and Ly, 

X,=X°, 

X,= X°, 

X,=x$, 

c=c°, 

P = P° 

After modifying equation (21) and using equation (19 and 20), the rate of 


change of ladius at any time t, in a control region is given by 



( 31 ) 


_ M S | C kc [0 ,X [a 
dt PsiC 

where, R a = Ro at time t = to and the subscript a represents the concentration and mole 
fraction of TMS in the /? u respectively 

With the radius in each f3 a at a given time, the pore volume (t) is 

computed and accordingly porosity is obtained by using equation (19) 

First the system of equation is solved forX,,X 2 ,andP, at each time step and 
then the iadius and the porosity of the pores at each control volume is computed using 
equation (31) This process is continued till the radius of the pore at the edges and 
corner of the domain under consideiation becomes zero 



CHAPTER 3 


RESULTS AND DISCUSSIONS 

The system of equations (27)-(29), in the preceding chapter, contained five 
variables X|, X 2 , X 3 , c and P The variable c and X 3 were eliminated from these using 
equations 24 and 30 

Thus, the general form of flux expression is, 

N=(A)(a l VX,+c/ 2 VX 2 +pVP), 

in which the coefficients (Xi, OC2 and P are dependent on the variables Xj, X2, X3, c and 
P First, we guessed the values for these variables to compute a’s and p and then 
substituted the flux expression, linear in the gradient of Xj, X 2 and P in equation (22), 
to obtain a system of finite difference equations A self consistent solution is 
iteratively obtained when the values guessed to compute a and p are the same as the 
solution of the finite difference linear equations 

Once the concentrations of all the gaseous species and pressure is computed 
everywhere at a given time, the deposition of SiC in a time step is computed by 
integrating equation 3 1 by the fourth order Runge-kutta method 

Repeated calculations of concentrations, and then marching forward in time 
determines the deposition profiles The process of densification ends when the radius 
of the pores at preform surface becomes zero 

The list of properties used in the calculations [34] are given below 



Table 1 List ot properties used in calculations 


Density of Silicon carbide, p, kg/m^ 

3210 

Activation energy, E A , kJ/mole 

120 

Pre-exponential factor, ko, m/s 

2 62 

Mole fraction of TMS, x° 

0 0452 

Mole fraction of HCL, x° 

0 

Pressure, P u , Pa 

1 0133 x 10 5 

! 


However, the first step in calculations is determining the anisotropy tensor, 
(A) , for all the volume elements 

3 1 DETERMINATION OF ANISOTROPY TENSOR 

As a sample calculation, consider a distribution of pores shown in Fig 3 1, 
where the bold lines define various (3 a and the lighter lines are the pores, the radius of 
the pores is the same in a but it is allowed to vary from one element to another 



Fig 3 1 Preform with pores in four orientations and divided in smaller 


elements (fla), Lj=L 2 =L 










Consider the / 3 a marked in the figure, where all pores lie in x-y plane Thus, 

» 

for all pores 6=90° and the pores are along <j> = 0°, 45°, 90°, and 135° with respect to 
the x-axis If we define ‘n’ to be the number of pores, with = 45° or 135°, below the 
longest pore in /? a (n=3 in Fig 3 1), then the total volume of the pores is 
V a p =2 (l + V2)(n + 1) 7tR 2 L, where R the radius of pores 


Hence the ODF for this pore arrangement is, 



and then using this ODF in the equation (14) of chapter 2, the anisotropy tensor 
obtained is, 



where s is the porosity of the element 

Similarly, (a) is computed for all the elements, and although, in time, the 

ODF does not change, a change in porosity forces re-computing (A) in every time 


step 


3 2 A TYPICAL CALCULATION 


The preform dimension in all cases is 1cm x 1cm The third dimension, used 
only for calculation of porosity, is arbitrary and taken as 0 001 cm, ensuring that it is 
larger than the initial pore radii, Ro=2 5 (im For the cases in which the orientations of 
the pores were symmetrical both with respect to the x- and y-axes, the calculations are 
only in the first quadrant of the preform Accordingly, the preform is descretized into 
51x51 node points and, when full preform is considered due to lack of symmetry, the 
number of node points is 101x101 

Various computer subroutines like NAG, IMSL were used to solve the system 
of equation The dnect subroutine using LU decomposition method, used from NAG 
consumed a lot of memory to solve the system of equations and, also, it took a long 
time The iteiative subioutines were found to be faster Finally a Bi-congugate 
giadient method was used to solve the linear system of equations However, the 
convergence was severely dependent on the initial guessed values To overcome this 
problem, a direct method was used only at the start of computations to provide good 
guesses for the subsequent iterations A typical calculation with 5 lx 51 mesh size 
lequired 11-12 Hours in Pentium-2 machine with 192MB RAM 

3 3 VALIDITY OF THE NUMERICAL APPROACH 

In older to evaluate the accuiacy of oui calculations, we devised two tests In 
the first, we take a simple pore structure, all pores are parallel to the x-axis, along (j) = 
0° This structuie is the first among the total of 1 1 orientations considered in this 
woik (see Appendix D) Foi this airangement of pores, the evolution of instantaneous 
porosity of the preform normalized by the initial porosity is shown in Fig 3 2 (01) 
However, since the pores are non-intersecting, and all of equal length, the results 


should be the same whether we apply this orientation averaging method in 
calculation, or consider just an isolated pore Therefore, to validate our calculations, 
we wrote a completely different, albeit much simpler, computer program, which takes 
the original dusty gas flux expression (without any orientation averaging), 
computes mass transport and deposition in a single pore The curve Ol in Fig 3 2 
was reproduced identically in this computation also 

Second, we also compared our results with those by Deepak and Evans [44] for 
a specific case, wheie pores are aligned along both x- and y-axes (see 06 in Appendix 
D) The noimali/ul poiosity as a function of piocessing time computed by the model 
given heie is computed with the results in [44] (see Fig 3 3) Clearly, the calculations 
are only slightly different at all times for all three temperatures The origin of 
difference lies m the manner the pore volumes have been taken 

Deepak and Evans consider only a specific pore structure, m which there are 
500 pores in both x- and y- directions in the first quadrant Because the pore structure 
if specified a prion, they are able to accurately include the volume of intersection of 
perpendicular pores and the total volume of the pores was 7 73x10 5 cm 3 Whereas, in 
the present work, since the model is for an arbitrary pore structure, we over-estimate 
the pore volume by counting volume of overlapping pores twice Thus, the initial 
pore volume for the same structure is 9 82x10 5 cm 3 , greater than the actual volume 
In fact, this disciepancy will be even larger when the number of pore intersections is 
greater Howevei, since a the enor in computing initial volume is also introduced in 
calculating pore volume subsequent time steps, the representation of normalized 
porosity, that is, e(t)/s 0 , is more accurate, accordingly, only a minor difference 
remains between the two calculations compared in Fig 3 3 







Foi this specific structure, if L, is the length of two intersecting pores 
associated with each intersection, then the ratio of actual porosity to that in this work 

Thus, the error between the two approaches is less if the pore radius 

is small or the length associated with each intersection is large (that is, less number of 
intersections), both of these imply lower volume of intersection Nonetheless, for 
pore radius 2 5 |im and starting porosity of the preform approximately 0 4, as that 
considered in this work, the error in calculations due to overestimating volume of 
intersecting pores is minimal 

3 4 EFFECT OF PORE ORIENTATION ON INFILTRATION 

In this work we have considered pores in four different orientations, 0 = 0°, 
45°, 90°, and 135° The various combinations include pores only in one, two, three or 
all four directions In all, there are a total of 11 pore networks, as presented in 
Appendix E, where, in addition to the pore network, we have given the ODF and the 
initial porosity The expressions for ODF, and corresponding porosity, are derived in 
terms on an integer ‘n’ i elated to the number of pores The integer ‘n’ is always odd 
and when an even number of pores is required, we use n+1 pores A value of ‘n’ is 
selected which gives the initial porosity closest to 0 4 Thus, we find the initial 
porosity of the preform considered in this work varies between 0 33 and 0 45, as 
represented in Table 2 
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Table 2 A comparison of infiltration behaviour of pores in various orientations 


Orientation 

Initial porosity 

Final porosity 

Final normalized 

number 

(eo) 

(e) 

porosity (e/e o) 

01 

0 393 

0 133 

0 337 

02 

0 389 

0 122 

0313 

03 

0 389 

0 100 

0 257 

04 

0 379 

0 117 

0 295 

05 

0 379 

0 102 

0 269 

06 

0 393 

0 114 

0 291 

07 

0 333 

0 082 


08 

0 402 

0 118 

HHESEBHI 

09 

0 402 

0 110 

0 274 

010 

0 451 

0 116 

0 257 

Oil 

0 379 

0 108 

0 284 


3 41 PARALLEL PORES 

First, we consider only the non-intersecting pores There are three such cases, 
01-3 In the first case, the pores are only along the x-axis, and symmetry requires 
calculation only in the first quadrant Because the prefoim we consider is square in 
cross-section, an equivalent structuie is obtained with all pores aligned along the y- 
axis In another case, the pores are along <j) = 45° However, when symmetry is 
assumed about both the axes, the structure (02) is different from the case when 
calculations are on full preform (03) The latter structure is also equivalent to <J) = 
135° 

We compare the progiess of infiltration for all three orientations in Fig 3 2 
Clearly, the pores onented 45° to the x-axis lead to greater densification 
Furthermore, a comparison of structures 02 and 03 intuitively suggest that the flux 































directed to the center of preform in 02 will lead to a greater reduction of porosity than 
in 03 In Figs 3 4 and 3 5, we have compared the normalized radii for the two cases 
at the end of the process, and similarly, in Figs 3 6 and 3 7, the normalized local 
porosity, which is equivalent to ratio of square of local radius divided by square of 
initial radius of the pores We find, contrary to expectation, the local radii and 
porosity, and final porosity are all lower for orientation 03 For 03, in quadrant 2 
and 4, this can be explained on the basis that flux is received in the pore from both the 
ends, as opposed to, only from one end in 02 However, when we compare Figs 3 8 
and 3 9 with Figs 3 4 and 3 6, respectively, we find even in quadiant 1 and 3, where 
the poie structure is identical in the two cases, the infiltration is greater for 03 But, 
in 02, while the porosity is uniformly distributed in all quadrants, m 03, porosity is 
lower in quadrants 2 and 4 compared to that in quadrants 1 and 3 And, in all 
quadrants, porosity is lowei for 03 than for 02 

Finally, densification is greater when the pores are inclined to the preform 
surface, lather than perpendicular to it (compare Ol to 02 and 03) 

3 4 2 PORES IN TWO ORIENTATIONS 

Fust, we have taken pores along (j) = 0°and 45° Assuming symmetry about 
both the axes, this structure is represented in 04, and if the structure is rotated by 90°, 
it is equivalent to poies along (j) = 45° and 90° Also, as in the case with single 
onentation pores, we have removed the symmetry requirement, for pores shown in 
05 In addition, we have considered pores that are either perpendicular (06) or 
inclined to the axes (07) In both cases, the structure is symmetric about the two 


axes 








y-axis 



Fig. 3.8: Only the first quadrant of Fig. 3.5 for comparison with Fig. 3.4. 
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When evolution of porosity during infiltration is compared in Fig. 3.10, there 
is practically no difference in results for 04 and 06. Furthermore, lowest porosity is 
shown for 07 

We also compare the local radii (Fig. 3.1 1-3.14) and porosity (Figs. 3.15-3.18) 
for the four orientations at the end of the process. For the same reason that causes 
greater infiltration m 03 compared to 02, a greater deposition occurs for 05 
compared to 04, as illustrated in Figs. 3.11 and 3 12 and Figs. 3.15 and 3.16. Again, 
even in the first quadrant, radii and porosity are lower for 05 (see Figs. 3.19 and 
3 20) However, the distribution of porosity in both 04 and 05 is not uniform across 
the area of the preform, whereas, in 06 and 07, the porosity contours are circular. 

Also, greater infiltration occurs when both pores are inclined rather than 
perpendicular to the preform surface (compare Fig. 3.13 with 3.17 and Fig. 3.14 with 
3.18). This, reinforces the fact that inclined pores lead to greater densification In 
other words, if fiber arrangement is such that it leads to perpendicular pores, then it 
may be prudent to carry out infiltration with pore structure as in 07. Then, after 
infiltration, rotate the preform by 45° and obtain the composite as in 06 by removing 
the extra triangular portions. This, of course, will take away the advantage of CVI, 
where near net shape composites can be processed. 

Finally, comparing porosity results for 01 and 06, we find for the same initial 
porosity, final porosity is lower when half the pores are in one orientation and other 
half perpendicular to it, rather than all pores in a single orientation 
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Time, Hours 


Fig 3 10: Evolution of porosity during densification of preform with 
pores in two orientations. 
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y-axis 



x-axis, cm 


Fig. 3.1 1: Final normalized radius for two orientation pores -- 04 
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x-axis, cm 

Fig. 3.17: Final normalized porosity for two orientation pores - 06. 
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Fig. 3-19: Final normalized radius for 05 only the first quadrant. 
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3.4.3 PORES IN THREE ORIENTATIONS 

In addition to two onentations parallel to the axes, vve also include pores 
oriented 45 to the x-axis With symmetry about both axes, the resulting structure is 
08; without the requirement of symmetry, the structure is as in 09. From Figs 3.21 
and 3 22, which are plots for local radii and, Figs. 3 23 and 3 24. for local porosity in 
08 and 09, respectively, we again find greater infiltration for 09: the final 
normalized porosity of the composite for 08 is 0.293, whereas it is 0 274 for 09. 

Also, as before, although the structure in the first quadrant is the same for the two 
cases, the local radii and porosity are lower for 09 (compare Fig 3 21 with 3.25 and 
Fig 3.23 with 3.26). 

Greater infiltration can be achieved by having more pores that are at an angle 
to the preform surface, rather than perpendicular to it. So, in Ol 0, out of three, only 
one set of pores are parallel to the x-axis and other two are in orientations 45°and 
135°. This structure is equivalent to one when pores are oriented as $ = 45°, 90°, and 
135°. According to results for local radii and porosity, m Figs. 3.27 and 3 28. 
respectively, the densification is greater than either 08 or 09. This fact is better 
illustrated in Fig. 3.29, where we show the evolution of porosity during the process. 

3.4.4 PORES IN FOUR ORIENTATIONS 

Finally, we present results (Figs. 3.30 and 3.31) for pores in all four 
orientations (Oil). Here, for a similar initial porosity as in 07, a part of porosity is 
redistributed along the two axes. Therefore, the infiltration is not as good as in 08. 
In fact, we have compared the final normalized porosity for all cases (see Table 2) and 
find that 08 yields largest densification. 
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In summary, the infiltration with single orientation and pores perpendicular to 
the preform surface is poorest; whereas most densification is achieved when only both 
the inclined pores are present In between are the cases, where a part or all the initial 
porosity due to inclined pores is redistributed between one or both the perpendicular 
pores 
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Fig. 3.24: Final normalized porosity for 09 
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Fig. 3.25. First quadrant radii in 09 for comparison with Fig. 3.21 
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Fig 3.27: Final normalized radius tor 010 
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CHAPTER 4 


CONCLUSIONS 


In the present work, a model for diffusion in porous media is developed. The 
porous media consisting of pore network oriented in specific directions are divided 
into a number of elemental volumes and by orientation averaging on the pore network 
in each elemental volume, the porous media is represented as a piece wise 
homogeneous medium through an anisotropy tensor The flux are based on Dusty gas 
model for describing mass transport in the porous medium. 

The present model is applied to study the effect of pore orientation on 
deposition of SiC in SiC preform by chemical vapor infiltration. The dusty gas based 
flux expression is used to generate a system of partial differential equations, 
containing variables mole fraction of TMS(Xi), mole fraction of HC1(X : ) and total 
pressure(P). The set of partial differential equations (PDE’s) with appropriate 
boundary conditions are solved, using finite difference method. The equations are 
discretized in a uniform grid and are solved by bi-congugate gradient method. 

The studies are carried out for pores oriented in four different orientations, <j) = 
0°, 45°, 90° and 135°. The various combinations included pores only in one, two, 
three and all the four directions. Total of 1 1 cases of pore orientations is considered. 
An initial pore radius of 2.5 (im is taken and the simulations are carried out at a 
temperature of 1200 K. 

For the pores oriented in only one direction, it is observed that the pores 
oriented at an angle of 45° through out the preform gives lowest porosity. Since for 
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this case, flux enters the pores from both the directions, as the pores are open in both 
the ends. 

For the case when the pores are oriented in two directions, the preform is 
densified more when the pores are inclined at angles of 0 = 45° and 135° to the 
preform surface This clearly shows that, if the pores are open m both the end for the 
flux to enter through them and also if they are inclined, then more densification of the 
preform can be achieved. Comparing the deposition for the case when the pores are 
oriented in single direction and in two directions, it is seen that more densification can 
be achieved when the pores are oriented as in the latter case. Also uniform deposition 
is observed when the pores are oriented perpendicular to each other (as in the case 
when pores are oriented at angles of <j) = 45° and 135° or <j) = 0° and 90°). 

For the pores oriented in three directions and all four directions. The same 
trend was seen as in the previous cases, that is, densification of the preform was 
dependent on the population of the pores inclined to the surface of the preform and 
also on the number of open ended pores. 

In brief it is seen that pores oriented perpendicular to the preform surface and 
in only one direction lead to higher porosity, whereas maximum densification 
occurred when the pores were inclined to the preform surface and oriented in two 
directions. 

Based on our calculation, we suggest that preform should be placed such that 
there are more number of open ended and inclined pores for the optimum 
densification of the preform. Then the densified preform can be cut to the desired 
shapes for commercial purpose. 
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half-thickness of ceramic preform in x and y direction, 
length, width and height of the elemental volume 
dyad used in the equation 
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APPENDIX A 



Figure A 1 The coordinates systems and pore orientation. 
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Fig. A.2: Schematic showing conversion from local to global coordinate system 
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APPENDIX B 


The anisotropy tensor due to orientation of the pores, which accounts for 

anisotropy of the diffusion phenomenon due to the orientation as well as the geometry 
of the pore is derived as, 


2717 : t 

(a) =— J J j PfcLM)/(e, <j>, ^)sin e dAie d<t> 

^000 

where the term /(0, §,£) is orientation distribution function 
The term M T LM is. 


^sin 2 0cos 2 <}> 
sin 2 0 cos <}) sin 
sin 0 cos 0 cos 0 

and 


sin 2 0 sin (J) cos <}) 
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For averaging at the region of control the delta functions were used, thus after 
averaging the anisotropy vector is calculated to be, 


f , , L, 

sin 0cos 0— 

£ 


2kk l 


(a)=— Iff 

w 47C J J J 
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APPENDIX C 

The Knudsen diffusion coefficient for species i is given by 



where m, is the molecular weight of gas species i 

The formula for the viscosity of a gas species i, T|„ is derived from kinetic 
theory of gases[43]* 

ri 1 =2.6693xl0' 5 

where p, is in gm cm' 1 s' 1 , m, is in gm, T in K and the lennard-jones parameter a, is 

in A 0 . Q 22 is the collision integral tabulated as a function of dimensionless 
k T 

temperature, T* =—2— , where e, is another Lennard-Jones parameter. Values of a, 

£ , 

and e, for H 2 and HC1 are known; but for TMS are not available in the literature. 
Therefore, they are estimated from emperical relationships[34], 

— =0 77T 
k B 

<7=0.841 V t ^ 

where the critical volume V c is in cm % lole » the critical temperature , T c is in K and 
computed and <7 are in K and A, respectively. 

Values of <7, and e, used in the computations are presented in Table B-l 


a, 2 a 22 (T‘) 
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Table B-T Lennard-Jones parameters for IMS, HC1 and H 2 


a, h 2 , a 0 ! 

2.915 

e/k B , H 2 , K 

38 

0 , HC1, A 0 

3.305 

e/k B , HC1, K 

360 

0 , TMS, A 0 

5.913 

£/k b , TMS, K 

398 


The viscosity of the gas mistime is given by Wilke's semi-emperical formula 


, - 1 


1=1 


where (p 






1 + 


^ J j 


m. 




f m vK 

1+B. 

V J 


The expression for the binary diffusion coefficient is also based on the kinetic 
theory of gases 


< 1_ + _P 


T 3 


1-3 


m, 


»»= L8583xl ° P 0 ;n"(T, : , 


where i3„ is in cm ~/ , T in K, m„ in gm, P in atm, 0 , 


f +0j ^ 


V “ J 


0 , k B T 

m A 0 and T„ t=— 5— 


where £„=^/e^ 
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APPENDIX D 


ODF AND ANISOTROPY TENSOR 


ORIENTATION 1 (01) 


y 



Fig. D.l: Schematic diagram of alignment of pore network in 01 


Fig. D 1 represents pore network for the case, w’hen the pores are aligned 
perpendicular with respect to z-axis (9 = 90°) and at an angle of 0 with respect to x- 

axis (<j> = 0°). The preform is divided into different elemental volume , where i 

can be 1 , 2, 3 or 4. The number 1 , 2, 3, and 4 represents the elemental volume at the 
interior, comers, lower and upper edges, and side edges of the preform respectively. 
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The dimension of the preform in the z-direction is taken as W for all the 
orientations. 

Anisotropy tensor for /? * 


Volume of the elemental volume is V a =LxLxW and total volume of the 


pores inside the elemental volume is. =jtR i (n + 1)L, where (n+1) represents 
number of pores in the element 

The ODF for the above pore network in the internal elemental volume is. 


above pore network is, 


VV 

V 2 J 


(8(<|))[(n + l)£8(t? - L)]} Thus, the anisotropy tensor for the 


(V 


^ 1 0^1 , (n + l)7tR- 

where e=2 - 

v° 


Anisotropy tensor for p 


LxW 

2 

a 


Volume of the element: V„ = — x— xW 

“ 2 2 


Total volume of the pores inside the elemental volume =7tR~(m + 1)— , for 


n - 1 

m= 

2 


The ODF: /(0,<M)= 


47t 

(m + l) L 
2 


7t 


8(<j))(m + l)£8 


e- 



The anisotropy tensor (a} = 



v 


0 ^ 


where e = 


(m + 1 )kR : 

/2 xW 


(n + l)7tR 2 
LxW 
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Anisotropy tensor for /? ^ 


Volume of element V =Lx— xW 

2 


Total volume of the pores: =7tR 2 (m + 1)L , where m = 


n-1 


The ODF is: f (&,§,£)= 


471 K ^ 


e- 


V 2 J 


{5((j>)[(m + \)i§{t - L)] } 


The anisotropy tensor. (a)=£ 


(m+l)L 

^1 (m + l)7tR : (n + l)7tR : 

, where, £ = —L = 4 ' 


v° °y 


Anisotropy tensor for /? 


a 


Ly/xW LxW 


Volume of the element: V„ =Lx— xW 

“ 2 


Total volume of the pores: V^ =7tR 2 (m + 1) L , where m ; 


n-1 


The ODF /(9,<M)= , 4 V §| 9~ 2(m + l)*8| 


( 71 ^ 


(m + l)L 


V 2 J 


(l 0) 


v 2 / 


v° °y 


The anisotropy tensor- (a) = s 

(m + 1 )tcR 2 (n + l)7iR 2 


Where, £ = 


L^xW LxW 
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ORIENTATION 2 (02) 



Fig. D.2: Schematic diagram of alignment of pore network in 02, (b = 45°. 

Anisotropy tensor for /?* 

Volume of the element is V a =LxLxW and total volume of the pores inside 
the elemental volume is V£ =V2 TtR 2 (n + 1)L , where n represents the number of 
pores oriented at an angle of (j) = 45° below the diagonal and is always odd. 

The ODF: 




-v/2 i 

M) 
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The anisotropy tensor (a)= e 


(2n + l) (2n + l) j 


2V2(n+l) 

(2n + 1) 


■: 1 , n-flK/27tR- 

, ■ . - , where, e = L 

(2n — 1 ) LxW 


2^2 


2V2(n + l) 2V2(n + l) j 


Anisotropy tensor for /? 


a 


Volume of the element: V a =— x— xW 

2 2 


Total volume of the pores inside the elemental volume: V£ =7tR 2 V2(m + 1)— , for 


n-1 

m = 

2 

The ODF 




471 


V2(m+ 1) — 
2 


7t ^ 

0-- l 

/. n ) 

to 

1 

1 

Ir 

m 

+ 2 Yt8 

1 

2 J 

r 4 ; 

L 

l 2 J 

jLmd 

1 = 1 

y lm+lJ2JJ 


( (2m + 1) 


The anisotropy tensor: (A)= £ 


2^2 {m -hi) 
(2m + 1) 

v 2V2(m + l) 


(2m + 1) ^ 

2V2(m + l) 

(2m + 1) 

2^2 (m + l) 


where, 


(m+l)V27iR 2 

L^xW 


(n +l)V2 7tR 2 
LxW 


Anisotropy tensor for p* 


Volume of the element: V a =Lx^xW 

2 

Total volume of the pores inside the elemental volume: V£=7tR 2 V2(m + l)L , where 


m 


n-1 

~2~ 
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The ODF: 


jTT^hr 5 e -v N o-t < m+2 ) <s F 2 2«s i-£^- 


■V2(im-1)L V 2A y 4 


2 “ i ' (m+1) 


' (3m +2) (3m +2) 

anisotropy tensor - ( a\= e 2-72 (m + l) 2^2 (m+ 1 ) j 
V ( 3m+2 ) (3m +2) 

4V2(m + l) 4^2 (m+l) I 


Where , e = 


(m+l)V2 7tR 2 _ (n+iyJlTcR 2 


LxW 


Anisotropy tensor for / 3 ‘ 


Volume of the element V„ =Lx — xW 


Total volume of the pores inside the elemental volume: V„ =7tR 2 V2(m + 1)L 


where, m 


The ODF: 


/(w)= 75^ 6 ri 


7t ^ / \ <■*( i — L ^ _ A 2, ,c r V2i L 
— (m+2)^6 £-42 — +2y\£& £—, -r— 

4/ ' t 2) tt 1 (m+1) 2 JJ 


f (3m +2) (3m +2) N 

T , . /aX 4V2(m + l) 4V2(m + l) 

The anisotropy tensor (a)= £ ( 3n !; +2 ) ( 3m +2) 

2-v/2(m + l) 2-v/2(m + l) 


Where, e 


(m + 1 yfl tcR ‘ (n + 1 yjl tcR ' 


LxW 
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k> | r 


ORIENTATION 3 (03) 


y 



Fig D.3 Schematic diagram of alignment of pore network 03 with (j> = 45° 
in entire preform. 

Fig. D.3 represents pore network similar to the 02, except symmetry 
requirement about both the axes is relaxed. 

The anisotropy tensor obtained for different elemental volumes are same as in 

02. 
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ORIENTATION 4 (04) 


y 



Fig. D.4 Schematic diagram of alignment of pore network in 04 <j ) = 0° 
and 45°. 

Anisotropy tensor for ^ 

Volume of the element is V a =LxLxW and total volume of the pores in the 
elemental volume is =(l+V2) (n + l)7tR' L , where n represents the number of 
pores oriented at an angle of <J> = 45°, below the diagonal pore and n+1 is the total 
number of pores parallel to the x-axis. 

The ODF: 


/(e,<M)= 


4xc 


IW2j(n+l)L 


x < 


K 


6(<J))[(n + - L)]+ §1 <J> - ^ 


re(<-^L)+2j<6| 


1=1 


vii T ^' 


(n+l) 
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The anisotropy tensor- (a)=£ 


(4n+3) 


(2n + l) ^ 


2(l + ^2](n-M) 2(1 + V2 )f n+l) i 

(2n+l) (2n-l) j 

2[l~+V2)(n-rl) 2(l + V2j(n+lJ j 


Where, 8- 


(n +l)(l + ^)nR 2 


LxW 


Anisotropy tensor for ^ 


Volume of the element: V„ =— x— x W 

2 2 


Total volume of the pores inside the elemental volume: V ! J='3tR 2 (l+V2)(m + l) 


where m 


The ODF: 


/(e,<M)= 


47t 

(l + V2)(mrH)| 


f 


0 -- 

2 


\ 

/ 


S(<|>) (m + l)f5 


L 'i 


^ 

2 


+ 6 




7C 


f8 


/ _ l a 




1 J 


+ 2^£8| 

1=1 


' V2i L ^ 

(m + l)2 




The anisotropy tensor. 




(A)« 




(4m +3) 

2(l + V2)(m + l) 
(2m +l) 

2(1 + V 2 |(m +l) 


(2m + 1) ' 

2(l + V2)(m + l) 
(2m +1) 

2(1 + V 2 )(m + 1) ^ 


Where, e= 



-xW 

2 


(n -t- 1 )(l + )tcR 2 
LxW 


Anisotropy tensor for 


Volume of the element: V 0 =Lx— xW 
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Total volume of the pores inside the elemental \olume: V==7tR : (l+-v / 2jfm + l)L, 
n -1 


where m = 


The ODF: 


/(0,<M)= 7 — 4 

(l + V2 j(m+l)L 


( ^ 7t ^ 
V 2 J 


8(ci))[(m + \)P.§{l - L)]+ 8 


f Tt ' 

0 -- 
V 4 A 


L „ v 2 ’ J , v/2 i L ^ 


(m + 2)£8' 1-42— + 2^£8 

1 2 1 


1=1 


/- 


(m + 1) 2 


JJ 


The anisotropy tensor • (a)=£ 


f (5m +4) (3m +2) ^ 

2(l + V2)(m+l) 2(1 + V2)(m + l) 


(3m +2) 


(3m +2) 




4(l + V2)(m + l) 4(l + V2)(m + l) 


Where, e= 


_(m+l)(l + V2)rtR 2 _(n+l)(l + V2)7tR 2 


L 


xW 


LxW 


Anisotropy tensor for /? 


a 


Volume of the elemental: V„ =Lx— xW 


Total volume of the pores inside the elemental volume: = 7 tR 2 (l+-\/ 2 )(m + l)L, 


where m = 


n -1 


The ODF: 


/(e,<M)= 


4ti 


(1 + V2j(m4- 1)L 


^0--^ 
V 2 2 


m 


2(m + 1)^8 


V 2 J. 


+ 8 


f 7t ' 

^-4 

V 4 4 


■L\ 




V2i L 


2-J 
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The anisotropy tensor (a}=£ 


f (7m +6) 

4(l + V2)(m ->-l) 

(3m +2) 

^2(1 + V2)(m + l) 


(3m -2 j 

4(1 + \/2 j(m + l) J 
(3m +2) i 

2|l+V2)(m + l) j 


where, c _ (^ + 0^ + )ttR 2 _ (n +l)(l + S)itR 


-xW 

2 


LxW 


ORIENTATION 5 (05) 


y 



Fig. D.5. 05 pore orientation; pores as in 04 without symmetry required 
both the axes. 

The anisotropy tensor obtained for different elemental volumes are 
same as in the previous case. 
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ORIENTATION 6 (06) 


y 



Fig. D.6: Schematic diagram of alignment of pore network in 06. <j) = 0° 
and 90°. 

Anisotropy tensor for * 

Volume of the element is V a =LxLxW and total volume of the pores inside 

the elemental volume is V„ = TcR ; 2(n + 1)L , where (n+1) represents number of pores 

in each direction. 

The ODF: 

/ (6,. 4>, eh 5^e -||8(<t.)[(n + l)<S(< - L)] + ^ - 1 j(n + !)«(<- L)]j . 
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The anisotropy tensor- (a)=s 


f\ \ 

- 0 ! 


0 - 
V 2 


, where, £ = 


2fn + lj7tR : 

Lx w 


Anisotropy tensor for ( 3 ^ 


Volume of the element: V =-x-xW 

2 2 


Total volume of the pores inside the elemental volume: = 7 tR 2 (m + l)L, where 


m = 


n-1 


The ODF: 


(m + 1 JL 


Vi' 

V 2 J 


8(<j) i (m + 1)£5| 


V 2 / 


+ 8 


f TC ^ 

^ 2 

V z A 


(m + l)f 8 


V 2 M 


The anisotropy tensor: (a)=£ 


A o' 

2 

0 I 

V 2 J 


where, £ = 


2(m + l)7tR 2 2(n + 1 )tiR 2 


%xW 


LxW 


Anisotropy tensor for /? 


a 


Volume of the element. V a =Lx— xW 


Total volume of the pores inside the elemental volume: V£ =7tR 2 2(m + 1) L , where. 


m =- 


n-1 


The ODF: 


/(e,<M)= 


471 


2(m + l)L 


C re ^ 

V 2 J 


( IT t 


8(<j))[(m + 1)^8(^ - L)]+ 8 


V 2 A 


2(m + l)£8j 


f t N 


V 2 J 
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T 0 1 

The anisotropy tensor. (a)=£ 2 ! Where, £ = 2fiTi-lj7iR 

o T I -X\V 


LxW 


Anisotropy tensor for /?' 


volume of the element: V =LX— xW 

2 


Total volume of the pores inside the elemental volume: =TtR 2 2(m + l)L, where. 


The ODF- 


/(e,^^)= 2 ( m 4 ; - 1 y^ 9 -||5W 2(m + l)^-|j +^-|j[(m + l^-L)]| 


/ \ ? u wu 2(m + l)7tR 2 2(n + l)7tR 

The anisotropy tensor. (a)=£ 1 , Where, £ = — — = — 

0 - -xW LxW 

2 2 
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ORIENTATION 7 (07) 


y 



D.7- Schematic diagram of alignment of pore network in 07. 0 = 45° and 135° 


Anisotropy tensor for p 


a 


Volume of the element is V a =LxLx W and total volume of the pores in the 
elemental volume is V a p = TtR 2 l4l (n + 1) L , where n represents the number of pores 
below the either diagonal. 


The ODF- 


/(e,<M)= 


4n 


2V2(n+ 1)L 


'a-*' 

v h 


x n 


A. 


C- 


+ 5 


<t>- 


3tc 


5(f-V2L)+2^f8 

1=1 

fS(f-V2L)+2^75 


a/ 2 i V 


( n + l) 


L 


v ("+» i 
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0 ! 
(2n-f-l) r 
2-\/2 ( n + 1) j 


Anisotropy tensor for ^ 

Volume of the element: V„ =— x— xW 

2 2 

Total volume of the pores inside the elemental volume: Y% = 7iR 2 2\/2(m + 1)^ , 


The anisotropy tensor: (a)=£ 


(2n + l) 
2-v/2(n+l) 

0 


where, g 


(n + 1)2V2 ttR 2 


LxW 


where, m = 


TheODF: 


/(e.<M)= 


471 


2V2(m+l)-= 1 


V + 


^ 7t ' 

V 4 T 


m 


1-^2 + 2£^5| 

V 1 J 


1=1 


l- 


+ 5 


f 


m 3% 

*"T 




f-4l- +2£f8j 
2 ; ■=> 


V2i L 

(m+l) 2 

, V2i L 

(m + l) 2 


The anisotropy tensor: (a)=e 


(2m + l) 
2^2 (m + l) 


(2m + 1 ) 
2^2 (m+l) 


where g 


(m+l)2V2rtR 2 (n +1)2^2 7tR 2 


L^xW 


LxW 


Anisotropy tensor for /P 


Volume of the element: V a =Lx— x W 
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Total volume of the pores inside the elemental \olume: V 0 ? =~R : 2\2fm + l)L, 
where, 


n -1 

m = 

2 

The ODF: 


/(e.<M)= 


471 


2V2(m+ 1)L 


7C 

2 j 


0 — - 


71 


V 


Jl 


+ sfo- — 

4 


l 


(m + 2x/ 1--I2- I+2V 

!, 2 J IT { (m + l)2 Jj 

/ r nr L ] rS r *\^2i L^j 

(m + 2K5| t - V2 - |+ 2X <5, < - ^ 


-7 


The anisotropy tensor. (a)=£ 


' (3m + 2) 

2^2 (m + l) 

0 


V 


0 




(3m + 2) 
4^2 (m + l) 


where e ( m + 1 W^R : _ (" + >W^ r! 

L ... LxW 

— xW 

2 


Anisotropy tensor for 


Volume of the element: V =Lx— x W 

2 

Total volume of the pores inside the elemental volume: =7tR 2 2V2(m + l)L, 


where, m 


n - 1 

1 
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The ODF. 


/(6,<M)= 


471 


f 7t ^ 


2V2(m+ 1)L 


V 2 y 


5 6- 


71 




~ L "i 


f, v'2i l' 1 


' 1 


+ *2X^£- h 

l 2 ) U { (m-rlj 2 


+ sf(j>-— (m-f-2)/5-G-V2— I+2Y /§i /- 7 — — 1 j 

l 4 j i 2j & ! (m + lj 2 JjJ 

^ (3m +2) 


The anisotropy tensor: (a) = e 


4^2 (m + l) 


0 


1 


(3m 4- 2) I’ 

2 V 2 (m + 1 ) j 


c _ (m + l)2V2 7tR 2 _ (n+l)2V27tR- 


-xW 


LxW 


ORIENTATION 8 (08) 



► x 


Fig. D.8: Schematic diagram of alignment of pore network in 08. <j> = 0° 
and 45°. 
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Anisotropy tensor for / ? 


a 


Volume of the element is V 0 — Lx Lx W and total volume of the pores m the 
elemental volume is V p = (2 + V2 ) (n + ljjtR 2 L, where n represents the number of 

pores oriented at an angle of <j) = 45°, below the diagonal. 

The ODF: 


4rt tO 


/(0, §,£)=-> 5| 

(2 + V2j(n+l)L 


v 2 y 


5(b)[(n + 1)/5(-4 _ l)]+ §f b - ^ [(n + 1)/5(/ - L)] 


1 


+ 5 


f . tz 
9 — 7 


-C 8 (jf — s/2 l) -i- 2 ^rn£ 


f V2 i 

! £-t— tL 





l V *JL 

1_ 


f (4n+3) 

(2n+l) ^ 

The anisotropy tensor: (a)=£ 

2(2 + -^2 )(n + 1 ) 
(2n+l) 

2(2 + V2)(n + l) 
(4n+3) 


v 2(2 + V2)(n + l) 

2(2 + V2)(n+l) y 


( n + l) jJJ 


where, g 


(n+l)(2 + V2)7iR : 

LxW 


Anisotropy tensor for /? 


a 


Volume of the element: V =— x— xW 

2 2 


Total volume of the pores inside the elemental volume: V a p =tcR 2 (2 + V^m + l)-^-, 


where m = — 
2 
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The ODF. 


/(e.<M)= 


4tc 


2 + V2)(m+ 1)- 


^ ( tc ^ 
V 2 y 


8(6 1 (m + 1 )>.8j -r 8 O- 

v 

-l n 


^ 1, 

(m-l/S £-— jji 

'v 2 tt\ 


( K N 


8 


m 


W 2 


The anisotropy tensor. (A)=e 


l v A v 

(4m +3) (2m -f-1) 

2(2 + V2)(m+l) 2(2 — V2j(m+1) 
(2m + l) (4m +3) 

2(2 + V2")(m + 1) 2(2 + V 2 j(m + 1) 
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Anisotropy tensor for 

Volume of the element V =Lx— xW 
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Total volume of the pores inside the elemental volume: V’=7tR 2 (2+72)fm + l)L, 
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ORIENTATION 9 (09) 


y 



Fig. D.9: Schematic diagram of alignment of pore network 09 with <j> = 0°, 
45° and 90° in entire preform. 

The anisotropy tensor obtained for different elemental volumes are 
same as in the previous case. 
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ORIENTATION 10 (010) 


y 



Fig. D.10: Schematic diagram of alignment of pore network in 010. (<j> =0°, 45°, 
and 135°. 

Anisotropy tensor for * 

Volume of the element is V a =LxLxW and total volume of the pores in the 
elemental volume is VP=(l+2V2)(n + l) 7 lR 2 L, where n represents the number of 
pores oriented at an angle of <j> = 45^ or 135°, below the either diagonal and (n+1 ) are 
pores parallel to x-axis. 
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where, £ = ( m + 1 )( 1 + 2 ^R 2 = (n ^ 2^2 IttR 2 

LxW ' 

Anisotropy tensor for f}^ 

Volume of the element: V =Lx-xW 
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Fig. D. 1 1 : Schematic diagram of alignment of pore network in 01 1 0 =0°. 45°, 90° 
and 135°. 


Anisotropy tensor for /?* 

Volume of the elemental volume is V a =LxLxW and total volume of the 
pores in the elemental volume is V„ =(2+2^2 ) (n + 1) 7tR" L , where n represents the 
number of pores oriented at an angle of (j) = 45° or 135°, below the either diagonal and 
(n+1 ) are the pores parallel to each axis. 

The 0DF: 
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Anisotropy tensor for ^ 


Volume of the element: V„ =Lx— xW 

“ 0 


Total volume of the pores inside the elemental volume: V p =7tR 2 (2+2-\/2)fm + l)L, 
n -1 


where, m 


The ODF- 


S(<J))[(m + 1)/5(/ - L)]+ S' 


V 1 A 


2(m + l)/S[ ] 


4tc 


/(9,<M)= 7 7 = r § 

2(1 + V 2 j(rn+ 1)L 


Q-*\U- n 


f -L > 


4 

3tc ^ 
V ^ ) 


{m + 2)ib \l-ji- +2]T£8 


+ 




J 


1=1 


(m + l) 2 


+ 


(m+2)fsff-V2- |+2^fS : 

l 2 J •=' 


\ V2i L N 

v (m + 1) 2 ^ 


The anisotropy tensor (a) = £ 


(4m + 3) 


3 


2(l + V2 )(m + 1 ) 

0 (5m+-4) 

4(1 + V2 )(m + 1 ) 


where, £ 


(m + 1)2(1 + V2)7tR : _ (n + 1)2(1 + V2)tiR 2 


-xW 

2 


LxW 


Anisotropy tensor for /?* 


Volume of the element - V =Lx— xW 

“ 7 


Total volume of the pores inside the elemental volume: 


V a p =7tR 2 (2+2V2)(m + l)L, 


where,, m = 


n - 1 

0 


105 



The ODF 


§(oj 2f m + 1 y o' Jf - — ,:fo o — — [i m - 1 • — L f]- 

i 1 ,J 




471 


2(1 + V2 )( 


m+ 1)L 


-i] 


5 


71 


V 4 t 


fm + 2)/6: /-V2— + 2Y /6 ■— 

' . 2 ) £ ; fm + lj 0 


S <j>- 


371 


fm + 2J/ s[ 7. - VI— |+ 2^ {Si / 


/-I ; 

' \2i L i] j 


) >=! I 


fm + lj 2 jj i 


The anisotropy tensor: (a)=£ 


(5m +4) 

4(l + V2)(m+l) 


0 


(4m +3) 

2(1 + -/2)(m + 1) 


where, £ = 


(m + l)2(l + V2)7tR 2 


— xW 
2 


(n +l)2(l + V2)7iR : 

LxW 


106 



LIST OF REFERENCES 


1 W . J. Lackey and T. L Starr, Fabrication of fiber reinforced ceramic composites 
by chemical vapor infiltration: Processing, Structure and properties.” pp. 39" -450 
m fiber reinforced composites. Edited by K. S Mazd'.yar.s: Noyes publications. 
Park Ridge, NJ, 1990. 

2. P. J. Lamicq, G. A Benhart, M. M. Dauchier, and J. G. Mace, “SiC SiC ccmncsitc 
ceramics,” Am. Ceram. Bull. Vol. 65, No 2, 1986, pp. 336-338 

3. E. Fitzer and R Gadow, “Fiber reinforced Silicon carbide.” Ceram Bull Vol 65. 
No 2, 1986, pp. 326-335. 

4 R Naslain, F Langlais, and R. Fedou, “The CVI -processing of ceramic matrix 
composities,” J De Physique, Colloque C5, supplement au n° 5, Tome 50, 1989. 
pp. c5-191-c5-207. 

5. G. S. Fishman and W. T. Petuskey, “Thermodynamic analysis and kinetic 
implications of chemical vapor deposition of SiC from Si-C-Cl-H gas system,” J 
Am. Ceram. Soc. Vol. 68, No. 4, 1985, pp. 185-190. 

6. Fitzer, W. Remmele, and S. Schoch, “Preparation of SiC-whiskers reinforced cvd- 
SiC,” J. De Physique, Colloque C5, supplement au n° 5, Tome 50, 1989, pp. c5- 
209-C5-217. 

7. R. D. Veltri, D. A. Condit, and F S. Galasso, “Chemical vapor deposited SiC 
matrix composites,” J. Am. Ceram. Soc. Vol. 72, No. 3, 1989, pp. 478-480. 

8. T. M Besmann, R. A. Lowden, D. P. Stinton, and T L. Starr, “A method for rapid 
chemical vapor infiltration of ceramic composites,” J. De Physique, Colloque C5, 
supplement au n° 5, Tome 50, 1989, pp. c5-229-c5-239. 


107 



9. D. P. Stmton, A. J. Caputo, and R. A Lou den, “Synthesis of fiber-reinforced SiC 
composites by chemical vapor infiltration,’' Am. Ceram. Bull. Vol. 66. No 2. 

1987, pp. 368-372. 

10 K Sugiyama and T Nakamura “Pulse CVI of porous carbon," J Mater Sci 
Letters. Vol 6, 1987, pp 331-333. 

1 1 M Sahimi, G. R. Gavalas, and T. T Tsotsis, “Statistical and continuum models of 
fluid-solid reactions m porous media,” Chem Engng. Sci Vol. 45, No 6,1990. 
pp. 1443-1502. 

12. S. Reyes and K. F. Jenson, “Estimation of effective coefficients m porous solids 
based on percolation concepts,” Chem Engng Sci Vol 40, No 9, 1985, pp 
1723-1734 

13. R Landaur, “The electrical resistance of binary metallic mixtures,” J App. Phvs 
Vol. 23, No. 7, 1952, pp.779-784. 

14. S. Kirckpatnck, “Classical transport in diordered media: Scaling and effectne 
medium theories,” Phys. Rev. Lett. Vol. 27, No. 25, 1971, pp. 1722-1725 

15. C. K. Harris, “Application of generalised effective-medium theory to transport m 
porous media,” Transport m porus media. Vol. , No. 5, 1990, pp. 517-542. 

16. V. N. Burganos and S. V Sotirchos, “Diffuison m pore networks: Effectn e 
medium theory and smooth field approximation,” AlChe Journal. Vol. 33, No. 10. 
1987, pp. 1678-1689. 

17. S. V. Sotnchos, “Multicomponent diffusion and convection in capillary structure," 
AlChe Journal. Vol. 35, No 12, 1989, pp. 1953-1961. 

18. P. G. Toledo, H. Ted Davis, and L. E. Scriven, “Transport properties of 
anisotropic porous media: Effective medium theory,” Chem. Engng. Sci. Vol 47. 
No. 2, 1992, pp. 391-405. 


108 



19 B. Evans, III, G M Watson, and E. A. Mason, “Gaseous diffusion m porous 
media at uniform pressure,”!. Chem. Phys Vol. 6, No. 12, 1961. pp 2076-20S3 

20. E A Mason, and A. P Malmauskas, “Gas transport in porous media, the dust\ 
gas model,” Chem. Eng. Monographs No 17. Else\ier Science Publishers, 1983 

21 R. Jackson, m Transport m porous catalysis, chem. Engng Monographs. Vol 4. 
1977. 

22 C. Feng and W. E Stewart, “Practical models for isothermal diffusion and flow of 
gases in porous solids,” Ind Eng. Chem. Fundam Vol 12.No. 2. 1973. pp. 143- 
147. 

23. E A. Efthimidiadis and S. V. Sotirchos, “A partially overlapping gram model for 
gas-solid reactions,” Chem. Engng. Sci. Vol 48, No 7, 1993, pp 1201-1212. 

24. J. Kim, J. Alberto ochoa, and S. Whitaker, “Diffusion in anisotropic porous 
media,” Transport m porous media. Vol , No. 2, 1987, pp. 327-356. 

25 M. P Hollewand and L. F Gladden, “Modelling of diffusion and reaction m 
porous catalyst using a random three-dimensional network model,” Chem. Engng 
Sci. Vol. 47, No. 7, 1992, pp. 1761-1770. 

26. M. P. Hollewand and L F. Gladden, “Representation of porous catalyst using 
random pore networks,” Chem. Engng. Sci. Vol. 47, No. 9-11, 1992, pp. 2757- 
2762. 

27. C. Rieckmann and F. J. Keil, “Multicomponent diffusion and reaction in three 
dimensional network. General kinetics,” Ind Eng. Chem Res. Vol. 36, 1997, pp. 
3275-3281. 

28. B. Fuertes and G. Marban, “Modelling gasification reactions including the 
percolation phenomenon,” Chem. Engng. Sci. Vol. 49, No. 22, 1994, pp. 3813- 
3821. 


109 



29 L. Zhang and N. A. Seaton, “The application of continuum eq_at:cr to diffusion 
and reaction m pore networks,” Chem. Engng. Sci Vol 49. Xo 1,1994. pp 41- 
50 

30 P. Rajniak and R. T. Yang, Linified network model for diffusion of condensable 
vapors in porous media,” AlChe Journal. Vol. 42. No 2. 1996. pp. 319-329 

31 S M. Gupte and J. A Tsamopoulos, “Densification of porous materials b\ 
chemical vapor infiltration,” J Electrochem. Soc. Vol 136, No 2, 1989, pp. 555- 
561. 

32 G. Q Lu, “Modelling the densification of porous structures m CVI ceramic 
composites processing,” Journal of Materials Processing Technolog} , Vol 37, 
1993, pp. 487-498. 

33. G. Y. Chung, B. J. McCoy, and M. Smith, “Chemical vapor infiltration: 
Modelling solid matrix deposition for ceramic composites reinforced with layered 
woven fabrics,” Chem Engng. Sci. Vol. 47, No. 2, 1992, pp. 311-323. 

34. D. Gupta and J. W. Evans, “Mathematical model for chemical vapor infiltration 
in a microwave-heated preform,” J. Am. Ceram. Soc Vol. 76, No 8, 1993, pp. 
1924-1929. 

35. M. Kawasa, Y. Ikuta, T. Tago, Takao Masuda, and Kenji Hashimoto, “Modelling 
of a thermal-gradient chemical vapor iniltration process for production of silicon 
carbide whisker/ Alumina composite,” Chem. Engng Sci. Vol. 49, No. 24A, 1994, 
pp. 4861-4870. 

36. H. Chang, T. F. Morse, and B. W. Sheldon, “Minimizing infiltration time during 
isothermal chemical vapor inilftration with methyltrichlorosilane,” J. Am. Ceram. 
Soc. Vol. 80, No. 7, 1997, pp. 1805-1811. 


110 



37. X. Tai, T Chou, and C M Ma, "‘Effects of deposition mechanisms m the 
modeling of Forced- flow.' /temperature —gradient chemical \ apor milftration." J 

Am. Ceram. Soc Vol. 77, No. 3, 1994, pp. 849-851 

i ! I ti | 

38 J Y Ofon and S. V Sotirchos, “Multiccmnc.'c.T mass transport m chemical 
vapor infiltration,” Ind Eng. Chem. Res. Vol. 35, 1996, pp 1275-1287 

39. J. Y Ofon and S V. Sotirchos, “Multidimensional modeling of chemical vapor 
infiltration, application to isobanc CVI,” Ind. Eng. Chem. Res Vol. 36, 1997, pp 
357-367. 

40. J. Y. Ofon and S. V. Sotirchos, “Optimal pressure and temperature for isobanc, 
isothermal chemical vapor infiltration,” AlChe Journal Vol 42. No 10, 1996, pp 
2828-2840. 

41. J. Y. Ofori and S. V. Sotirchos, “Structural model effects on the predication of 
chemical vapor infiltration models,” J. Electrochem Soc. Vol. 143, No. 6, 1996, 
pp. 1962-1973. 

42 J. O. Herschfelder, C. F. Curtiss and R. B Bird, m Molecular Theory' of Gases and 
Liquids, John Wiley, New York, 1964. 

43. R. B. Bird, W. E. Stewart, and E. N. Lightfoot, in Transport Phenomena, Wiley, 
New York, 1960. 

44. J. W. Evans and D. Gupta, “ Mathematical modelling of the production of 
composites by CVI,” in Materials Processing in computer Age, Ed. V. R Voller, 
M. S. Stachowicz and B. G. Thomas, Proceedings Of Int. Symp. Sponsered by 
TMS, Feb 17 -Feb -21, 1991, New Orleans, pp 111-122 


111 



1 311847 



a t30SS7 

P 

This book is to be returned on the 
date last stamped. 




